A Comparison between a Minijet Model and a Glasma Flux Tube Model for 
Central Au-Au Collisions at v /s/\T/v=200 GeV. 
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In is paper we compare two models with central Au-Au collisions at ^/sjvat=200 GeV. The first 
model is a minijet model which assumes that around ~50 minijets are produced in back-to-back pairs 
and have an altered fragmentation functions. It is also assumed that the fragments are transparent 
and escape the collision zone and are detected. The second model is a glasma flux tube model 
which leads to flux tubes on the surface of a radial expanding fireball driven by interacting flux 
tubes near the center of the fireball through plasma instabilities. This internal fireball becomes an 
opaque hydro fluid which pushes the surface flux tubes outward. Around ~12 surface flux tubes 
remain and fragment with ~l/2 the produced particles escaping the collision zone and are detected. 
Both models can reproduce two particle angular correlations in the different pti Pti bins. We also 
compare the two models for three additional effects: meson baryon ratios; the long range nearside 
correlation called the ridge; and the so-called mach cone effect when applied to three particle angular 
correlations. 



PACS numbers: 25.75.Nq, 11.30.Er, 25.75.Gz, 12.38.Mh 
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I. INTRODUCTION AND REVIEW OF 
MODELS 

In this paper we discuss two models. The first 
model is a minijet model [ll. The second is a glasma 
flux tube model (GFTM)[|. 

The paper is organized in the following manner: 
Sec. 1 is the introduction and review of models. 
Sec. 2 discuss two particle angular correlation in the 
two models. Sec. 3 discuss baryon and anti-baryon 
formation in both models. Sec. 4 demonstrates how 
the ridge is formed by flux tubes when a jet trigger 
is added to the GFTM. Sec. 5 treats the so-called 
mach cone effect by analyzing three particle angular 
correlations in the two models. Sec. 6 presents the 
summary and discussion. 



A. Minijet Model 

The analysis of angular correlations led to unan- 
ticipated structure in the final state of p-p and 
Au-Au collisions, subsequently identified with par- 
ton fragmentation in the form of minijets [3-7]. 
Two-component analysis of p-p and Au-Au spec- 
tra revealed a corresponding hard component, a 
minimum-bias fragment distribution associated with 
minijets, suggesting that jet phenomena extend 
down to 0.1 GeV/cjH From a given p-p or Au- 
Au collisions particles are produced with three kine- 
matic variables (rj,(p,pt)- The p t spectra provide 
information about parton fragmentation, but the 
fragmentation process is more accessible on a loga- 
rithmic momentum variable. The dominant hadron 
(pion) is ultrarelativistic at 1 GeV/c, and a rela- 



tivistic kinematic variable is a reasonable alterna- 
tive. Transverse rapidity in a longitudinally comov- 
ing frame near midrapidity 77=0 is denned by 



y t = ln([m t +p t ]/m ), (1) 

with mo a hadron mass and m 2 = p\ + m\ . If one 
integrates over <j) the event multiplicity {dn/ drj = p, 
the ID density on rj) can be written for p-p collision 
in a two component model as 



p{yf,v = 0) = So(yt) +H (y t ) 



(2) 



So is a Levy distribution on m t which represents 
soft processes and Hq is a Gaussian on y t which rep- 
resents hard processes. The soft process has no <f> 
dependence but does have a Gaussian correlation in 
the longitudinal direction^). This correlation can 
be expressed in a two particle way as Arj = 771 — 772 
which is the difference of the psuedorapidity. 



p s (A V ) = A exp[-Ar] 2 /2a 2 } 



(3) 



The hard process arise from a scattering of two 
partons thus minijets are formed. Each minijet frag- 
ments along its parton axis and generates a 2D cor- 
relation A<p = 4>x — 4>2 and A 77 = rji — 772. 



p h (A<t>,Ari) = A h exp[-Act> 2 /2o-l\exp[-Aif /2al]. 

Since for every minijet fragmenting leading to a 
peak at A<fi = 0°, there is its scattered partner in the 
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backward direction A<fi = 180°. The backward scat- 
tered minijet will range over many psuedorapidity 
values, therefore its correlation with the fragmenta- 
tion of the near side minijet will have a broad A-q 
width. In this condition simple momentum conser- 
vation is good enough and that is a —cos(A<p) term. 

The two-component model of hadron production 
in A-A collisions assumes that the soft component is 
proportional to the participant pair number (linear 
superposition of N-N collisions) , and the hard com- 
ponent is proportional to the number of N-N binary 
collisions (parton scattering) |l0| . Any deviations 
from the model are specific to A-A collisions and 
may reveal properties of an A-A medium. In terms 
of mean participant path length v = Inun I n par t the 
Pt-integrated A-A hadron density on 7/ is 



dn 



Tipart drjAA 



= Ps 



vph- 



(5) 



By analogy with Eg. (2) A-A density as a function 
of centrality parameter v becomes 



-Paa{uu v,r] = G) = S NN (y t )+vr AA (y t ; v)H NN {y t 



tpart 



(6) 

In the above equation S/viv(yt) = Sb(yt) an d 
HffN(yt) — Ho(yt) which is p-p scattering or 
nucleon-nucleon scattering. We also can define for 
the hard density for A-A collisions in terms of v 
as H AA (y t \v) = r AA (y t \v)H NN (y t ). The density 
2/n par t 1/27T l/yt d 2 /dy t dr] for pions as a function of 
Vt at V s aw = 200 GeV for p-p, Au-Au, and ifjvjv(yt) 
is shown in Figure 1. 

We see from the above equation the SVivCj/t) is 
universal and scales with the participant pairs. This 
means we can extract vH AA (yt) from the density 
measured in p-p and Au-Au at y / J/vw=200 GeV 
from Figure 1 giving us Figure 2. 

Finally the ratio H aa /Hmn (Raa) is plotted in 
Figure 3. In the Au-Au central collision (0-12%) at a 
yt value of 2 (p t =0.5 GeV/c) there is 5 times as many 
pions coming from minijet fragmentation as there is 
in a N-N collision. This implies a large increase in 
correlated pion fragments and should show up as an 
increased in two particle angular correlations. We 
will see this in Sec. 2 where we show these angular 
correlations and discuss the number of particles in 
the minijets. 

The common measurement of parton energy loss 
is expressed as a nuclear modification factor R AA . 
The hard-component ratio r AA measured over y t 
(pt) provides similar information. Note that y t = 
5 {pt— 10 GeV/c) r AA = 0.2 which is the same 
value calculated for R AA . This value is considered 



caused by jet quenching as partons are absorbed by 
the opaque medium. However in the minijet pic- 
ture suggests that no partons are lost in A-A col- 
lisions. Their manifestation (in spectrum structure 
and correlations) is simply redistributed within the 
fragment momentum distribution, and the fragment 
number increases. A high-p t triggered jet yield may 
be reduced by a factor of five within a particular p t 
range, but additional fragments emerge elsewhere, 
still with jet-like correlation structure [H 1(1 fllj. 



B. Glasma Flux Tube Model 

A glasma flux tube model (GFTM)0 that had 
been developed considers that the wavefunctions of 
the incoming projectiles, form sheets of color glass 
condensates (CGC)[l2| that at high energies collide, 
interact, and evolve into high intensity color elec- 
tric and magnetic fields. This collection of primor- 
dial fields is the Glasma 0, Q], and initially it is 
composed of only rapidity independent longitudinal 
color electric and magnetic fields. An essential fea- 
ture of the Glasma is that the fields are localized 
in the transverse space of the collision zone with a 
size of 1/Qs- Qs is the saturation momentum of par- 
tons in the nuclear wavefunction. These longitudinal 
color electric and magnetic fields generate topologi- 
cal Chern-Simons charge[l5j which becomes a source 
for particle production. 

The transverse space is filled with flux tubes of 
large longitudinal extent but small transverse size 
^Qj 1 . Particle production from a flux tube is a 
Poisson process, since the flux tube is a coherent 
state. The flux tubes at the center of the trans- 
verse plane interact with each other through plasma 
instabilities [H, [l6[ and create a locally thermalized 
system, where partons emitted from these flux tubes 
locally equilibrate. A hydro system with transverse 
flow builds causing a radially flowing blast wave [13] • 
The flux tubes that are near the surface of the fire- 
ball get the largest radial flow and are emitted from 
the surface. 

Q s is around 1 GeV/c thus the transverse size of 
the flux tube is about l/4fm. The flux tubes near 
the surface are initially at a radius ~5fm. The (f> 
angle wedge of the flux tube is ~l/20 radians or 
^3°. Thus the flux tube initially has a narrow range 
in </>. The width in Ai] correlation of particles results 
from the independent longitudinal color electric and 
magnetic fields that created the Glasma flux tubes. 
In this paper we relate particle production from the 
surface flux tube to a related model Parton Bubble 
Model(PBM)[H|. It was shown in Ref.pJ] that for 
central Au-Au collisions at % /ijv r /v=200 the PBM is a 
good approximation to the GFTM surface flux tube 
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FIG. 1: Pion yt density for five Au-Au centralities (solid curves). Density for p-p is also shown (solid dots). Hnn is 
the hard Gaussian on y t which represents hard processes (dash-dot curve). 



formation. 

The flux tubes on the surface turn out to be on the 
average 12 in number. They form an approximate 
ring about the center of the collision see Figure 4. 
The twelve tube ring creates the average behavior of 
tube survival near the surface of the expanding fire 
ball of the blast wave. The final state surface tubes 
that emit the final state particles at kinetic freezeout 
are given by the PBM. One should note that the 
blast wave surface is moving at its maximum velocity 
at freezeout (3c/4). 

The space momentum correlation of the blast 
wave provides us with a strong angular correlation 
signal. PYTHIA fragmentation functions [20J were 
used for the tube fragmentation that generate the fi- 
nal state particles emitted from the tube. The initial 
transverse size of a flux tube ~l/4fm has expanded 
to the size of ^2fm at kinetic freezeout. Many par- 
ticles that come from the surface of the fireball will 



have a pt greater than 0.8 Gev/c. The final state 
tube size and the Hanbury-Brown and Twiss (HBT) 
observations [2l| of pions that have a momentum 
range greater than 0.8 GeV/c are consistent both 
being ~2fm. A single parton using PYTHIA forms 
a jet with the parton having a fixed rj and <f> (see Fig- 
ure 5). For central events each of the twelve tubes 
have 3-4 partons per tube each at a fixed <f> for a 
given tube. The p t distribution of the partons is 
similar to pQCD but has a suppression at high p t 
like the data. The 3-4 partons in the tube which 
shower using PYTHIA all have a different rj values 
but all have the same </> (see Figure 6). The PBM 
explained the high precision Au-Au central (0-10%) 
collisions at ^/saw = 200 GeV[22j](the highest RHIC 
energy). 
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FIG. 2: The hard component of pion y t spectra in the form vHAA(yt) (thicker curves with changing line style) 
compared to two-component reference uHNN(yt) (dotted curve). The dashed reference curves are limiting cases for 
v = 1, 6. 



II. THE CORRELATION FUNCTION FOR 
CENTRAL AU-AU DATA 

We utilize a two particle correlation function in 
the two dimensional (2-D) space 1 of A<fr versus At). 
The 2-D total correlation function is defined as: 

C(Ad>, A V ) = S(A<j>, Ar])/M(Acj), A v ). (7) 

Where S(A<fi, Arj) is the number of pairs at the 
corresponding values of A<p, Arj coming from the 



1 A</> = </>i — <f>2 where </> is the azimuthal angle of a particle 
measured in a clockwise direction about the beam. A77 = 
r)i — r]2 which is the difference of the psuedorapidity of the 
pair of particles 



same event, after we have summed over all the 
events. M(A<j), A-q) is the number of pairs at the cor- 
responding values of A(f>, A-q coming from the mixed 
events, after we have summed over all our created 
mixed events. A mixed event pair has each of the 
two particles chosen from a different event. We make 
on the order of ten times the number of mixed events 
as real events. We rescale the number of pairs in the 
mixed events to be equal to the number of pairs in 
the real events. This procedure implies a binning 
in order to deal with finite statistics. The division 
by M(A0, At;) for experimental data essentially re- 
moves or drastically reduces acceptance and instru- 
mental effects. If the mixed pair distribution was the 
same as the real pair distribution C(A</>, Arj) should 
have unit value for all of the binned Acj), Arj. In 
the correlations used in this paper we select parti- 
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FIG. 3: Hard-component ratio for pions and five Au-Au centralities (thicker curves with changing line styles) relative 
to the N-N hard-component reference. The connected dots are data from NSD p-p collisions. 



cles independent of its charge. The correlation of 
this type is called a Charge Independent (CI) Cor- 
relation. This difference correlation function has the 
denned property that it only depends on the differ- 
ences of the azimuthal angle ( A<fi) and the beam an- 
gle (Ary) for the two particle pair. Thus the two di- 
mensional difference correlation distribution for each 
tube or minijet which is part of C(A0, Arj) is similar 
for each of the objects and will image on top of each 
other. We further divide the data (see Table I) into 
Pt ranges (bins). 



pt range 


amount 


4.0GeV/c - 


- l.lGeV/c 


149 


l.lGeV/c- 


- Q.8GeV/c 


171 


0.8GeV/c - 


OmGeV/c 


152 


0.65GeV/c 


- 0.5GeV/c 


230 


0.5GeV/c - 


- OAGeV/c 


208 


OAGeV/c - 


- 0.3GeV/c 


260 


0.3GeV/c - 


- Q.2GeV/c 


291 



TABLE I: The pt bins and the number of charged parti- 
cles per bin with \r)\ < 1.0. 



Since we are choosing particle pairs, we choose for 
the first particle pn which could be in one bin and 
for the second particle pti which could be in another 
bin. Therefore binning implies a matrix of pn vs pt%. 
We have have 7 bins thus there are 28 independent 



combinations. Each of the combinations will have a 
different number of enters. In order to take out this 
difference one uses multiplicity scaling [23, 24]. The 
diagonal bins one scales event average of Table I. For 



Plane section of tube geometry perpendicular to the beam at eta = 



G 




FIG. 4: The tube geometry is an 8 fm radius ring perpendicular to and centered on the beam axis. It is composed 
of twelve adjacent 2 fm radius circular tubes elongated along the beam direction as part of the flux tube geometry. 
We project on a plane section perpendicular to the beam axis. 



the off diagonal combinations one uses the product 
of square root of corresponding diagonal event aver- 
ages. In Figure 7 we show the correlation function 
equation 7 for the highest diagonal bin p t 4.0 to 1.1 
GeV/c. Figure 8 is the smallest diagonal bin p t 0.3 
to 0.2 GeV/c. 

Once we use multiplicity scaling we can compare 
all 28 combinations. In Figure 9 we show the 28 plots 
all having the same scale. This make it easy to see 
how fast the correlation signals drop off with lower- 
ing the momentum. These plots show the properties 
of parton fragmentation. P1P7 has the same signal 
as P2P6, P3P5, and P4P4. It should be noted that 
both the minijet model and GFTM give the same 
two particle correlations. 



variable amount fluctuations 



minijets 


48 


4 


particles 


13 


4 


soft 837 29 



TABLE II: Parameters of the minijet model for charged 
particles. 

tation (see Table II). All of the particles coming from 
minijet fragmentation add toward the final observed 
correlation signal with none being absorbed. The 
fact that the spectrum has been soften and spread 
out in the beam direction (77), is a medium modifica- 
tion which has not been yet calculated using QCD. 



A. The properties of the minijet model 

For the above correlations on average there was 48 
minijets per central Au-Au collision. Each minijet 
on the average showered into 13 charged particles. 
The soft uncorrelated particles accounted for 837 or 
57% of the charged particles. This means that 43% 
of the charged particles come from minijet fragmen- 



B. The properties of the Glasma Flux Tube 
Model 

For the above correlations on average there was 
12 final state tubes on the surface of the fireball per 
central Au-Au collision. Each tube on the average 



7 



Jet 



Parton Shower 




Parton 

Fixed (j) and rj 



FIG. 5: A jet parton shower. 



showered into 49 charged particles. The soft un- 
correlated particles accounted for 873 or 60% of the 
charged particles. Since the tubes are sitting on the 
surface of the fireball and being push outward by 
radial flow, not all particles emitted from the tube 
will escape. Approximately one half of the particles 
that are on the outward surface leave the fireball 
and the other half arc absorbed by the fireball (see 
Figure 10). This means that 20% of the charged par- 
ticles come from tube emission, and 294 particles are 
added to the soft particles increasing the number to 
1167 (see Table III). 

The particles that emitted outward are boosted in 
momentum, while the inward particles are absorbed 
by the fireball. Out of the initial 49 particles per 
tube the lower p t particles have larger losses. In Ta- 
ble IV we give a detailed account of these percentage 
losses and give the average number of charged par- 
ticles coming from each tube for each p t bin. 



variable 


amount fluctuations 


tubes 


12 


particles 


24.5 5 


soft 


1167 34 



TABLE III: Parameters of the GFTM for charged par- 
ticles. 



In the surface GFTM we have thermalization and 
hydro flow for the soft particles, while all the two 
particle angular correlations come from the tubes 
on the surface. The charge particle spectrum of the 
GFTM is given by a blastwave model and the di- 
rect tube fragmentation is only 20% of this spec- 
trum. The initial anisotropic azimuthal distribution 
of flux tubes is transported to the final state leaving 
its pattern on the ring of final state flux tubes on 
the surface. This final state anisotropic flow pattern 
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Fixed (j) Different T| 

FIG. 6: Each tube contains 3-4 partons as shown. 



Pi 



amount %survivc 



A.OGeV/c 


- l.lGeV/c 


4.2 


100 


l.lGeV/c 


- Q.8GeV/c 


3.8 


76 


Q.8GeV/c - 


- OmGeV/c 


3.2 


65 


OMGeV/c 


- 0.5GeV/c 


4.2 


54 


0.5GeV/c 


- QAGeV/c 


3.2 


43 


OAGeV/c 


- Q.3GeV/c 


3.3 


35 


0.3GeV/c 


- Q.2GeV/c 


2.6 


25 



TABLE IV: Parameters of the GFTM for p t of the 
charged particles. 



can be decomposed in a Fourier series fa , Vn, v$, 
...). These coefficients have been measure [25j and 
have been found to be important in central Au-Au 
collisions. We will come back to these terms later on 
when we consider three particle angular correlations. 



III. FORMATION OF BARYON AND 
ANTI-BARYON IN BOTH MODELS 



It was shown above that charged particle produc- 
tion differ between p-p and Au-Au at ^snn — 200. 
High pt charged particles in Au-Au are suppressed 
compared to p-p for p t > 2.0 GeV/c (y t > 3.33 for 
pions). Figure 3 shows taa for pions which has a 
suppression starting at y t — 3.0 (p t — 1.5 GeV/c). 
This shift between suppression of charged particle 
at pt > 2.0 GeV/c and suppression of pions at pt 
= 1.5 GcV/c is made up by an enhancement of 
barvonsfiZrjj. Proton taa (see Figure 12) which is 
the most numerous baryon shows an enhancement 
starting at p t > 1.5. This enhancement continues to 
Pt = 4.0 (y t = 4.0) 2 . Finally at p t =10.0 GeV/c (y t 



Note yt is calculated using a pion mass. 
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FIG. 7: A(j> vs. At; CI correlation for the 0-5% centrality bin for Au-Au collisions at ^snn 
particles to be in bin 7 pt greater than 1.1 GeV/c and pt less than 4.0 GeV/c. 



200 GeV requiring both 



5-0) t~aa = 0.2 which is the same value as pions. 



B. Glasma Flux Tube Model 



A. Minijet Model 



There is no a priori model for the proton hard 
component. The excess in the proton hard com- 
ponent appears anomalous, but may be simply ex- 
plained in terms of parton energy lossQ. 



We have shown previously that inside the tube 
there are three to four partons with differing longi- 
tudinal momenta all at the same <f>. The pt distri- 
bution of the partons inside the tube is similar to 
pQCD but has a suppression in the high p t region 
like the data. The partons which are mainly gluons 
shower into more gluons which then fragment into 
quarks and anti-quarks which overlap in space and 
time with other quarks and anti-quarks from other 
partons. This leads to an enhance possibility for 
pairs of quarks from two different fragmenting par- 
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FIG. 8: A<j) vs.Ar] CI correlation for the 0-5% centrality bin for Au-Au collisions at ^/sjvjv = 200 GeV requiring both 
particles to be in bin 1 p t greater than 0.2 GeV/cand p t less than 0.3 GeV/c. 



tons to form a di-quark, since the recombining par- 
tons are localized together in a small volume. The 
same process will happen for pairs of anti-quarks 
forming a di-anti-quark. This recombination process 
becomes an important possibility in the GFTM com- 
pared to regular jet fragmentation. Since the quarks 
which overlap have similar phase space, the momen- 
tum of the di-quark is approximately twice the mo- 
mentum of the quarks, but has approximately the 
same velocity. When mesons are formed quarks pick 
up anti-quarks with similar phase space from frag- 
menting gluons to form a color singlet state. Thus 
the meson has approximately twice the momentum 



of the quark and anti-quark of which it is made. 
When the di-quark picks up a quark and forms a 
color singlet it will have approximately three times 
the momentum of one of the three quarks it is made 
from. Thus we expect a pt spectrum scaling when 
we compare mesons to baryons. Figure 12 shows 
the ratio protons plus anti-protons to charged parti- 
cles as a function of pt for particles in our simulated 
central Au-Au collisions. In Figure 12 we also plot 
the ratio from central Au-Au RHIC data[26]. These 
experimental results agree well (considering the er- 
rors) with the GFTM predictions for all charged par- 
ticles. The background particles which came from 
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FIG. 9: A(j> vs. At/ CI correlation for the 0-5% centrality bin for Au-Au collisions at ^snn ~ 200 GeV for all 28 
combination of pti vs pt2 not rescaled (see text). 



HIJINGJI^I have the same ratios observed in p-p col- 
lisions, while particles coming from our tube have a 
much larger ratio. 



IV. THE RIDGE IS FORMED BY THE 
FLUX TUBES WHEN A JET TRIGGER IS 
ADDED TO THE GFTM 

In heavy ion collisions at RHIC there has been 
observed a phenomenon called the ridge which has 
many different explanations[28-35]. The ridge is a 
long range charged particle correlation in A77 (very 



flat), while the A<j) correlation is approximately jet- 
like (a narrow Gaussian). There also appears with 
the ridge a jet-like charged-particle-pair correlation 
which is symmetric in At? and A</> such that the peak 
of the jet-like correlation is at A77 = and A</> = 0. 
The A0 correlation of the jet and the ridge are ap- 
proximately the same and smoothly blend into each 
other. The ridge correlation is generated when one 
triggers on an intermediate pt range charged particle 
and then forms pairs between that trigger particle 
and each of all other intermediate charged particles 
with a smaller p t down to some lower limit. 
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In the rest frame of the bubble 1/2 
of the particles are emitted away 
from the fireball surface 



average of 49 particles 24.5 go outward 



FIG. 10: Since the tubes are sitting on the surface of the fireball and being push outward by radial flow, not all 
particles emitted from the tube will escape. Approximately one half of the particles that are on the outward surface 
leave the fireball and the other half are absorbed by the fireball. 



A. STAR experiment measurement of the 
ridge 

Triggered angular correlation data showing the 
ridge was presented at Quark Matter 2006 [36]. Fig- 
ure 14 shows the experimental A0 vs. A?y CI cor- 
relation for the 0-10% central Au-Au collisions at 
^/snn = 200 GeV requiring one trigger particle pt 
between 3.0 and 4.0 GeV/c and an associated par- 
ticle pt above 2.0 GeV/c. The yield is corrected for 
the finite A77 pair acceptance. 

In this paper we will investigate whether the 



GFTM can account for the ridge once we add a jet 
trigger to GFTM generator [? ]. However this trig- 
ger will also select jets which previously could be ne- 
glected because there was such strong quenching[37- 
39] of jets in central collisions when a jet trigger had 
not been used. We use HIJING|27j to determine the 
number of expected jets. We have already shown 
that our final state particles come from hadrons at 
or near the fireball surface. We reduce the num- 
ber of jets by 80% which corresponds to the esti- 
mate that only the parton interactions on or near the 
surface are not quenched away and thus at kinetic 
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FIG. 11: Hard-component ratio for protons and five Au-Au centralities (thicker curves with changing line styles) 
relative to the N-N hard-component reference. The features are only comparable to Figure 3 for 60-80% or y t > 4.5. 
Note y t is calculated using a pion mass. 



freezeout fragment into hadrons which enter the de- 
tector. This 80% reduction is consistent with sin- 
gle 7T° suppression observed in Ref . [391] . We find for 
the reduced HIJING jets that 4% of the Au-Au cen- 
tral events (0-5%) centrality at ^/snn = 200 have a 
charged particle with a pt between 3.0 and 4.0 GeV/c 
with at least one other charged particle with its p t 
greater than 2.0 GeV/c coming from the same jet. 
With the addition of the quenched jets to the gener- 
ator, we then form two-charged-particle correlations 
between one-charged-particle with a pt between 3.0 
to 4.0 GeV/c and another charged-particle whose pt 
is greater than 2.0 GeV/c. The results of these cor- 
relations are shown in Figure 15. 

We can compare the two figures, if we realize that 
the away-side ridge has around 420,000 pairs in Fig- 



ure 14 while in Figure 14 the away-side ridge has 
a correlation of around 0.995. If we multiply the 
correlation scale of Figure 14 by 422,111 in order 
to achieve the number of pairs seen in Figure 13, 
the away-side ridge would be at 420,000 and the 
peak would be at 465,000. This would make a good 
agreement between the two figures. We know in our 
Monte Carlo which particles come from tubes which 
would be particles form the ridge. The correlation 
formed by the ridge particles is generated almost en- 
tirely by particles emitted by the same tube. Thus 
we can predict the shape and the yield of the ridge 
for the above pt trigger selection and lower cut, by 
plotting only the correlation coming from pairs of 
particles that are emitted by the same tube (see Fig- 
ure 15). 
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FIG. 12: Shows the ratio of protons plus anti-protons to charged particles as a function of pt for particles in our 
simulated central Au-Au collisions. We also plot the ratio from central Au-Au RHIC data. These experimental 
results agree well (considering the errors) with the GFTM predictions for all charged particles. The plotted ratio for 
the background particles coming from HIJING is similar to p-p collisions. 



In Ref. [36| it was assumed that the ridge yield was 
flat across the acceptance while in Figure 14 we see 
that this is not the case. Therefore our ridge yield 
is 35% larger than estimated in Ref. [36| . Finally we 
can plot the jet yield that we had put into our Monte 
Carlo. The jet yield is plotted in Figure 16 where we 
subtracted the tubes and the background particles 
from Figure 14. 
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^/iivjv = 200 GeV requiring one trigger particle pt between 3.0 to 4.0 GeV/c and an associated particle pt above 2.0 
GeV/c. 



Triggered particles 




FIG. 14: The GFTM generated CI correlation for the 0-5% centrality bin requiring one trigger particle p t above 3.0 
GeV/c and less than 4.0 GeV/c and another particle p t above 2.0 GeV/c plotted as a two dimensional A(j> vs. Arj 
perspective plot. The trigger requirements on this figure are the same as those on the experimental data in Figure 
13. 
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Ridge 




FIG. 15: The ridge signal is the piece of the CI correlation for the 0-5% centrality of Figure 15 after removing all 
other particle pairs except the pairs that come from the same tube. It is plotted as a two dimensional A(j> vs. Arj 
perspective plot. 



Jet Signal 




FIG. 16: The jet signal is left in the CI correlation after the contributions from the background and all the bubble 
particles are removed from the 0-5% centrality (with trigger requirements) of Figure 15. It is plotted as a two 
dimensional vs. A77 perspective plot. 
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B. PHOBOS experiment measurement of the 
ridge 

The PHOBOS detector select triggered charged 
tracks with p t > 2.5GeV/c within an acceptance of 
< jf" 9 < 1.5. Associated charge particles that es- 
cape the beam pipe are detected in a range \rj\ < 3.0. 
The CI correlation is shown in Figure 17 [40]. The 
near-side structure is more closely examined by inte- 
grating over | Acj)\ < 1 and is plotted in Figure 18 as a 
function of A77. PYTHIA simulation for p-p events 
is also shown. Bands around the data points rep- 
resent the uncertainty from flow subtraction. The 
error on the ZYAM procedure is shown as a gray 
band at zero. All systematics uncertainties are 90% 
confidence level. 

We generate using our above GFTM the two- 
charged-particle correlations between one-charged- 
particle with a pt between greater than 2.5 GeV/c 
which has an < 0.75 and another charged-particle 
whose p t is greater than 7 MeV/c in a range |7y| < 
3.0. The results of this correlation is shown in Fig- 
ure 20. We see that the triggered correlation is vary 
similar to the PHOBOS results. In order to make 
a comparison we integrate the near-side structure 
over \A(f>\ < 1 and plot it as a smooth curve in Fig- 
ure 21 as a function of Arj. We also again plot the 
PHOBOS points on the same plot. The long range 
correlation over Arj produced by the GFTM is pos- 
sible and does not violate causality, since the glasma 
flux tubes are generated early in the collision. The 
radial flow which develops at a later time pushes the 
surface tubes outward in the same <f> direction be- 
cause the flow is purely radial. In order to achieve 
such an effect in minijet fragmentation one would 
have to have fragmentation moving faster than the 
speed of light. 



V. THE MACH CONE EFFECT AND 
THREE PARTICLE ANGULAR 
CORRELATIONS 

It was reported by the PHENIX Collaboration^ 
that the shape of the away-side Acj> distributions 
for non-peripheral collisions are not consistent with 
purely stochastic broadening of the peripheral Au- 
Au away-side. The broadening and possible change 
in shape of the away-side jet are suggestive of a mack 
cone (43. The mach cone effect depended on the 
trigger particles and the associated particles used 
in the two particle angular correlations. The mach 
cone shape is not present if one triggers on 5-10 
GeV/c particles and also use hard associated parti- 
cles greater than 3 GeV/c[43|. The away-side broad- 
ening in the two-particle depends on the v-i subtrac- 



tion and the method of smearing the away-side dijet 
component and momentum conservation. 

One needs to go beyond two particle correlations 
in order to learn more. A three particle azimuthal 
angle correlation should reveal very clear pattern 
for the mach cone compared to other two particle 
method. The STAR Collaboration Q has made 
such correlation studies and does find a structure 
some what like a mach cone. However the trigger de- 
pendence as seen in Ref.[43[ and the measurement 
of the conical emission angle to be independent of 
the associated particle pt in Ref.[H| is not consis- 
tent with a mach cone. The similarity of the mach 
cone and the ridge is very interesting to consider and 
makes one consider that they are the same effect. In 
the last section we showed that flux tubes can ex- 
plain the ridge, thus they should explain the mach 
cone. One would not expect the minijet model will 
give a mach cone like correlation. 

We saw at the end of Sec. 2 that final state 
anisotropic flow pattern can be decomposed in a 
Fourier series {v\, i>2, W3, ■•■)• When one triggers 
on a particle coming from a flux tube, the other flux 
tubes contribute the components of the anisotropic 
flow pattern. In Figure 22 we show two such tube ar- 
rangements. In order to make contact with the data 
and show the difference between the minijet model 
and GFTM we will define a trigger particle and a 
reference particle(s). We choose a p t of greater than 
1.1 GeV/c for the trigger and less than this value for 
the reference. The two particle correlation for this 
trigger and reference for the minijet model and the 
GFTM are equal to each other. 

We define a three particle correlation using the 
azimuthal angle between a trigger particle 1 and a 
reference particle 2 (A0i2) and for the third particle 
the azimuthal angle between a trigger particle 1 and 
a reference particle ?>(A4>i^). The two particle corre- 
lation of the two models are the same and the three 
particle correlation is nearly the same. In order to 
obtain the true three particle effect we must remove 
the two particle correlation from the raw three par- 
ticle correlation. This removal gives the so-called 
three particle cumulant. Figure 22 and Figure 23 
shows the three particle cumulant for A0i2 vs. A^>i3 
for the minijet model and the GFTM. 

The minijet model Figure 22 shows only diagonal 
away-side response coming from the underlying dijet 
nature of the minijets. The GFTM Figure 23 also 
shows the diagonal away-side response of momentum 
conservation, however ^3 configurations of Figure 21 
gives an off-diagonal island at A</f)i 2 ~ 130°, A0i 3 ~ 
230° or vice versa. 

Let us look at the reported STAR data[45j in or- 
der to address this off-diagonal effect. In Figure 24 
we show the raw three particle results from STAR 
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FIG. 17: PHOBOS triggered CI correlation from 0-30% central Au-Au collisions. It is plotted as a two dimensional 
A<f> vs. A77 perspective plot. 



where (a) is the raw two particle correlation (points) . 
They also show in (a) the background formed from 
mixed events with flow modulation added- in (solid) . 
The background subtracted two particle correlation 
is shown as an inset in (a) where the double hump 
mach cone effect is clear. In Figure 24(b) the STAR 
raw three particle correlation is shown which re- 
quired one trigger particle 3.0 < p t < 4.0 GeV/c 
and other reference particles 2 and 3 with 1.0 < pt 
< 2.0 GeV/c. An off-diagonal island does appear at 
A0i2 ~ 2.3, A0i3 ~ 4.0 radians or vice versa. This 
is the same values as the island which occurs in the 
GFTM and is the off-diagonal excess that is claimed 
to be the mach cone. Like the ridge effect the mach 
cone seems to be the left overs of the initial state flux 
tube arrangements related to the fluctuations in the 
third harmonic. Higher order harmonic fluctuations 
become less likely to survive to the final state. 



VI. SUMMARY AND DISCUSSION 

In this article we have made a comparison between 
two very different models for central Au-Au colli- 
sions. Both models are successful at describing the 
the spectrum of particles produced and the two par- 
ticle angular correlations observed in ultrarelativis- 
tic heavy ion collisions. The first model is a minijet 
model which assumes that around ^50 minijets are 



produced in back-to-back pairs and has an altered 
fragmentation function from that of vacuum frag- 
mentation. It also assumes that the fragments are 
transparent and escape the collision zone and are 
then detected. The second model is a glasma flux 
tube model which leads to longitudinal color electric 
and magnetic fields confined in flux tubes on the sur- 
face of a radial expanding fireball driven by interact- 
ing flux tubes near the center of the fireball through 
plasma instabilities. This internal fireball becomes 
an opaque hydro fluid which pushes the surface flux 
tubes outward. Around ~12 surface flux tubes re- 
main and fragment with ~l/2 the produced particles 
escaping the collision zone and are detected. 

We expand our comparison to other phenomenon 
of the central collisions. We considered in Sec. 
3 baryon and anti-baryon formation in both mod- 
els. There was no a priori reason for the ex- 
cess in the minijet model, while in the glasma flux 
tube model(GFTM) recombination of quarks into di- 
quarks and anti-quark into anti-di-quarks leads to a 
natural excess of baryon and anti-baryon formation 
in this model. 

The formation of the ridge phenomena is discussed 
in Sec. 4. In order to achieve a long range cor- 
relation effect in minijet fragmentation one would 
have to have fragmentation moving faster than the 
speed of light. The GFTM however can have a long 
range correlation over Arj since the glasma flux tubes 
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_ Near-side, \Aty\ < 1 .0 




FIG. 18: PHOBOS near-side trigger yield integrated over | Arj}\ < 1.0 for 0-30% central Au-Au compared to PYTHIA 
p-p (dashed line) as a function of Arj. Bands around the data points represent the uncertainty from flow subtraction. 
The error on the ZYAM procedure is shown as a gray band at zero. All systematics uncertainties are 90% confidence 
level. 



are generated early in the collision. The radial flow 
which develops at a later time pushes the surface 
tubes outward in the same <fi direction because the 
flow is purely radial. Thus a long range A77 last to 
the final state of the collision. A very good compar- 
ison was achieved between data and GFTM. 

Sec. 5 treats the so-called mach cone effect by 
analyzing three particle angular correlations in the 
two models. The minijet model and the GFTM have 
the same two particle angular correlations but when 
the three particle azimuthal angular correlations are 
compared the two models differ. The minijet model 
shows only diagonal away-side response coming from 
the underlying dijet nature of the minijets, while 
GFTM also shows a diagonal away-side response it 
also shows an off-diagonal island. Like the ridge ef- 
fect the mach cone seems to be left over from the ini- 
tial state flux tube arrangements related to the fluc- 
tuations in the third harmonic(i>3). This off-diagonal 
island excess is seen in the data and is claimed to be 
the mach cone. 

Relativistic Heavy Ion Collider (RHIC) collisions 



are conventionally described in terms of two major 
themes: hydrodynamic evolution of a thermal bulk 
medium and energy loss of energetic partons in that 
medium through gluon bremsstrahlung. The minijet 
model is not consistent with is standard view. The 
glasma flux tube model generates a fireball which 
becomes an opaque hydro fluid that is consistent 
with conventional ideas. Even though both mod- 
els can obtain the same spectrum of particles and 
the same two particle angular correlations, it is only 
the GFTM that can tie all together. 
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These events add to Mach Cone 

..... 




These events add to v 2 
These are substracted away 



FIG. 21: Top arrangement has a V3 component while the bottom gives a 172 component. In most mach cone analyses 
V2 is subtracted. 
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minijet cumulant 




FIG. 22: The minijet generated three particle cumulant for the 0-5% centrality bin requiring one trigger particle 
Pt above 1.1 GeV/c and other reference particles p t below 1.1 GeV/c plotted as a two dimensional combinations of 
trigger particle 1 paired with two reference particles 2 and 3 creating A<^>i2 vs. A</>i3 perspective plot. 
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flux tube cumulant 




FIG. 23: The GFTM generated three particle cumulant for the 0-5% centrality bin requiring one trigger particle 
Pt above 1.1 GeV/c and other reference particles p t below 1.1 GeV/c plotted as a two dimensional combinations of 
trigger particle 1 paired with two reference particles 2 and 3 creating A<^>i2 vs. A</>i3 perspective plot. 
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FIG. 24: (a) STAR data raw two particle correlation points, background from mixed events with flow modulation 
added-in (solid) along with the background subtracted two particle correlation (inset), (b) Raw three particle corre- 
lation from STAR requiring one trigger particle 3.0 < pt < 4.0 GeV/c and other reference particles 2 and 3 with 1.0 
< Pt < 2.0 GeV/c creating a two dimensional perspective plot A(/>i2 vs. A0i3. 



